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I. INTRODUCTION 

The sun is a source of copious fluxes of 
charged particles which escape into inter- 
planetary space. These particles range in 
energy from the few keY solar wind parti- 
cles to the several hundreds of MeV parti- 
cles produced by the larger solar flares. 
There is even growing evidence that the 
sun may be a nearly continual emitter of 
low energy (several tens of keV) protons. 
This review is limited in that it concerns 
itself essentially entirely with the cha- 
racteristics of the solar particles accel- 
erated by solar flares and subsequently 
observed near the orbit of the earth. 

The number of solar flares and the 
fluxes of energetic 20 MeV) solar 
particles observed at the earth varies 
in a manner similar to that of the sunspot 
number during the eleven-year solar cycle. 
This is Illustrated by the data of Pig. 1 
where the smoothed sunspot numbers for 
cycles 19 and 20 are plotted as a function 
of time. Also shown are histograms of the 
yearly integrated proton intensities for 
protons £ 30 MeV for both solar cycles 
(A. J. Masley, private communication). 
These particle fluxes are obtained from 
riometer measurements of solar proton- 
produced PCA events in the polar-cap re- 
gions. The solar particle fluxes peaked 
in total intensity a year or more after 
the sunspot maximum during cycle 19* It 
remains to be seen if this same phenomenon 
holds during the current cycle. 

This review discusses in order solar 
particle intensity-time profiles, the com- 
position and spectra of solar flare events, 
and the propagation of solar particles in 
interplanetary space. The last section. 


dealing with the effects of solar parti- 
cles at the earth, discusses riometer ob- 
servations of polar cap cosmic noise ab- 
sorption events and the production of so- 
lar cell damage at synchronous altitudes 
by solar protons . 

IX. INTENSITY- TIME PROFILES 


Detectability Limits 

The first observations of energetic par- 
ticles due to solar production were the 
sea-level measurements of Lange a nd For- 
bush (1942), and Forbush ( 1946) . Using^ 
shielded ionization chambers built to ob- 
serve galactic cosmic rays, large enhance- 
ments in the chamber counting rates on 2o 
February 1942 (~1 day prior to a large 
magnetic storm), on 7 March 1942, and 25 
August 1946 were observed. The develop- 
ment of the cosmic ray neutron monitor in 
the late 1940 f s and the super neutron mo- 
nitor in the late 1950 ' s and early I960' s 
have enabled many more solar particle in- 
creases to be observed on the ground. The 
neutron multiplicity monitor (Nobles et 
al., 1967) enables particle spectral in- 
formation to be obtained from a single 
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Other pre-spacecraft observations of 

olar cosmic rays were made during the IG 
erlod by Polar cap radio absorption tech- 
iques ( Bailev . 1957 ) and by balloon mea- 
urements"! The balloon observations of 
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— — t-ttt a. j -p ool q rirntnnS . 


INTEGRATED SOLAR PROTON INTENSITY DISTRIBUTION 
_ SUNSPOT CYCLES 19 AND 20 

tr 



193 



The neutron monitor has continued to be 
a valuable tool for the detection of ener- 
getic solar particles. The world-wide 
deployment of stations provide data for 
studying both the direction of incidence 
of the primary particles as well as their 
energy. An example of the difference in 
response to an event by super neutron mo- 
nitors at two different energy (or rigidi- 
ty) cut-off latitudes is shown in Fig. 2 
for data measured during the 28 January 
1967 solar event (Bukata et al., 1969 ) . 

The vertical cut-off rigidity for the 
Churchill station is 1.0 GV (determined 
essentially entirely by the atmospheric 
cut-off) while that of Dallas, at mid- 
latitudes geomagnetically, is 4.35 GV. 



Fig. 2 Neutron monitor observations at Ft. Churchill and 
Dallas of the 28 January 1967 solar event 


The time-intensity profile of the 
Churchill monitor response to the January 
1967 event is quite typical of the classi- 
cal, diffusive-like profiles recorded by 
high energy flare-particle detectors and 
will be discussed in more detail in Sec- 
tion V. In general in diffusive-like 
events, a rapid rise to the peak particle 
intensity is followed by a slower, ex- 
ponential or power- law decay with time. 

With the advent of instrumentation flown 
on spacecraft, the energy sensitivity 
threshold for the detection of solar parti- 
cles was dramatically reduced. This reduc- 
tion in the lower limit 0 f the energy of 
particle detectability has continued until 
today measurements of 300-500 keV solar 
protons are routinely carried out. Accom- 
panying this decrease in the energy sensi- 
tivity of particles that can be measured 
was an increase in the number of solar 
events that were observed. Furthermore, 
at the lower energies measured, the time 
histories of the events became complex 
with no simple relationships often evident 
between events. 


The data plotted in Fig. 3 illustrates 
the enormous differences in the descrip- 
tion of the interplanetary particle inten- 
sities that could be made during a one- 
month period depending upon the energy 
sensitivity limits available for analysis. 
The data were obtained by the solar proton 
monitoring experiment on the Explorer 4l 
satellite (C. 0 . Bostrom, private communi- 
cation). If only the higher energy channel 
(E > 60 MeV) were available for analysis, 
only one large event (March 27) and a small 
event (March 24) would have been apparent. 
Although the decay times are longer, both 
of these events had a diffusive temporal 

appearance similar to the neutron monitor 
event of Fig. 2. 

As the particle energy threshold in Fig. 

3 is lowered, more solar events are detect- 
ed. The event beginning on March 6 (when 
viewed in the E > 10 MeV channel) no longer 

has a diffusive shape. The solar fluxes 
in the 1-10 MeV channel of Fig. 3 are ob- 
served to remain above their background 
level through out the entire 31 -day period 
plotted. 



MARCH , 1970 , UT 


Figure 3. Solar protons measured in inter- 
planetary space by the solar proton 
monitoring experiment on Explorer 41 
during March 1970. Data courtesy 
of C. 0. Bostrom. 
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Not all of the interplanetary particle 
enhancements result from discrete flare 
events. In general, flare-associated 
events occur in close association with so- 
lar X-ray and microwave emissions. Fur- 
ther, as noted above, the time-intensity 
profiles of flare-associated events tend 
to have a diffusive appearance. Three 
other types of particle enhancements, in 
addition to the flare-associated events 
considered in this paper, have been classi- 
fied and discussed extensively in the li- 
terature. These are: 

a) Particles associated with active 

centers: The onsets of these par- 

ticles at the earth display no velo- 
city dispersion and appear to be 
co-rotating with solar-active cen- 
ters. Such enhancements have been 
observed to occur each solar rota- 
tion for many successive rotations 
( e.g. , Fan et al. , 1968; McDonald 
and Desai , 1971 ) • 

b) Recurrent events: These particle 
increases occasionally occur in the 
next so'lar potation following a 
flare. They appear to originate 
from the same active region as that 
producing the flare ( e.g. , Bryant 
et al . , 1965)- 

c) Energetic storm particles: Enhance- 

ments of low energy protons that 
appear for several hours around the 
time of occurrence of interplanetary 
shock waves ( e.g. , Axford and Reid , 
1963 ; Bryant et a 1.7 1965; Rao et 

al. , 1967)7 Proton enhancements 

lasting for several minutes, appar- 
ently resulting from acceleration 
at the shock front, have been re- 
ported ( e.g. . Singer , 1970; Lanzer^ 
otti, 19"69a7 1970a; Armstrong and 
Krimigis, 1970; Ogllvle and Arens , 
1971)7“ 

Although the three solar particle en- 
hancements listed above are important for 
understanding solar processes and inter- 
planetary propagation, they will not be 
elaborated upon here. 

Data Organization 

Although the intensity- time profiles of 
high energy flare particles are similar in 
their overall diffusive appearance, abso- 
lute differences as a function of particle 
energy are common. Cline and McDonald 
( 1968 ) have shown that the time history 
of the high energy proton and electron 
fluxes from the 7 July 1966 solar flare 
is dependent upon particle velocity. This 
is evident in Fig. 4 where the observed- 
time profiles for three proton and one 
electron channel are plotted in Fig. 4a. 

In Fig. 4b, the four particle flux channels 
have been normalized to their peak values 
and the abscissas have been transformed to 
represent the distance traveled from the 
flare occurrence. 


59 TO 80 MeV PROTONS 
38 TO 59 McV PROTONS 
16 TO 38 MeV PROTONS 


<R =■ 330 TO 400 MV , 0 - .34 TO .39) 
(R ^ 260 TO 330 MV , fi = .28 TO .34) 
(R = 175 TO 260 MV, ^ .18 TO .28) 



(a) 


(b) 


Fig. U (a) Profiles of the observed intensity of protons (16-38 MeV, 

38-59 MeV, and 59-81 MeV) and electrons {>3 MeV) plotted versus 
tine following the T July 1966 flare; (b) profiles of each particle 
channel relative intensity (i/l ) plotted as a function of distance 
traveled ( x *vt) where v is the B mean velocity for each particle 
channel ( Cline and McDonald , 1968) . 


Although the higher energy particle 
fluxes from this flare can be organized 
quite well by considerations of velocity- 
dependent travel, Lin (1970a) has shown 
that when electrons of energy > 45 keV 
from this event are included in the anal- 
ysis, they show a broader curve than those 
in Fig. 4b. Furthermore, the E > 45 keV 
electrons appear to arrive earlier than 
the protons and electrons considered by 
Cline and McDonald ( 1968 ). Lin and Ander - 
son (1967) and Lin ( 1970a) have interpreted 
this earlier arrival to low energy electron 
production either higher in the solar at- 
mosphere or prior to the proton production. 

III. COMPOSITION 

The most recent reviews of soiar cosmic 
ray composition are those of Biswas and 
Fichtel (1965) and Fichtel (1970)7 They 
discuss in detail the several counter and 
emulsion measurements made on balloons and 
rockets beginning during the maximum of 
solar cycle 19. The discussion here will 
be limited primarily to observations of 
solar alpha particles and electrons. Ob- 
servations of higher-Z elements will be 
briefly outlined. 
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Solar Alpha Particles 


The primary characteristic arising from 
the observations of solar alpha particles 
is that the ratio of the fluxes of solar 
alphas to solar protons appears to vary 
widely between individual events and even 
within a single event. Both of these char- 
acteristics can be seen from Fig. 5 (Durga- 
prasad et al M 1967). Here are plotted the 
proton to alpha ratios for several differ- 
ent events as a function of particle kinet- 
ic energy. 

More recently, using satellite instru- 
mentation, the solar alpha measurements 
have been extended to lower energies and 
the time resolution during a single event 
has been substantially improved (Armstrong 
et al. , 1969; Lanzerotti and Robbins . 197 b). 
A comparison of the intensity- time profiles 
of solar protons and alphas from the series 
of flares on 21 and 23 May 1967 , is shown 
in Fig. 6 ( Lanzerotti and Robbins , 1970 ). 
The overall appearance of the intensity 
profiles of the two species are similar 
although differences do exist. In parti- 
cular, the energetic storm particle en- 
hancement at the time of the sudden com- 
mencement (SC) on May 24 is not strongly 
evident in the alpha fluxes. 

The detailed alpha to proton ratios 
throughout the May 23 event are shown in 
Fig. 7 . Large changes in the ratios are 
observed, particularly at the low energies, 
for protons and alphas when compared as to 
equal energy and equal energy per nucleon 




Pig. 6 Alpha particle and proton fluxes measured following 
the 21 and 23 May 1967 solar events. ( Lanzerotti 
and Robbins. 1969). 



KINETIC ENERGY ( MEV/ NUC) 

Fig. 5 Ratio of protons to hell* nuclei as a function of kinetic energy per nucleon 
measured at different tines in several solar events. The events in which the 
measurements were made are: A, B -12 November i 960 ( Biswas et. al . . 1962 ); 

C, D, E - 15 November 196o( Biswas et. al . , 1963 ); F - 3 ~ 5gpwifre r i 9 60 ' ' 

( Biswas and Fichtel, 1964); G - 2 6 March 1964 and 5 February 1965 (McDonald et. 

_ , -.t Wo5(; P - 2 September 1966 ( Dnrgaprased et. al.. 1967). (From ' 

EAirgaprased et . al . . 196?). 
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PIr. 7 Alpha particle to proton flu* ratios 

(equal velocities) end equal particle 
energy per charge (Lanzerotti and 
Robbins, 1969). 
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(equal velocity). However, the central 
panel of Fig. 7 indicates that after the 
increases in the ratios following the May 
24 sudden commencement, the ratios remain 
constant throughout the remainder of the 
event for particles compared as to equal 
energy per charge. Lanzerottl and Robbins 
( 1970 ) have interpreted the increases in 
the alpha to proton ratios observed after 
the sudden commencement on May 24 as a 
source effect. They suggest that the so- 
lar particles observed prior to the SC 
were predominantly from flares on May 21 
whereas those particles observed after the 
May 24 increase were from flares on May 23 . 
They also suggested that the constancy of 
the ratios for equal energy per charge may 
indicate an important role for electric 
fields in low energy particle propagation 
and/or acceleration. Similar behavior of 
the alpha to proton ratios following other 
flare events have been noted ( Lanzerottl , 
1970a; Lanzerottl and Graedel ; 1970)^ 

Heavy Nuclei 

Heavy solar cosmic ray nuclei (Z>3) were 
first detected in nuclear emulsion stacks 
flown on a rocket during the 30 September 
i 960 event ( Flchtel and Guss , 1961) . The 
evidence gained from a number of balloon, 
rocket, and satellite experiments in the 
early I 960 ’ s indicates that the spectral 
forms for solar heavy nuclei and solar 
alpha particles are the same in any one 
event for particles down to ~30 MeV/nuclecn 
(Flchtel, 1970). 

A statistical study using satellite data 
has been made of the ratio of solar alphas 
to Z>3 nuclei for a number of events in 
1967 - 1968 . The study indicates that for 
particles of E > 0.5 MeV/nucleon, the 
spectral behavior observed at the higher 
energies continues to hold ( Armstrong and 
Krimigis , 1971). It was found that the 
event-integrated alpha to heavy ratio was 
^20 ± 10 for most events, a value substan- 
tially smaller than the ratio of 48 ± 8 
reported by Durgaprasad et al . ( 1968 ) after 
the 2 September 1966 event in the 12-35 
MeV/nucleon range. It Is also smaller 
than the weighted mean of 58 ± 5 deter- 
mined from six large events in 196 O-I 969 
( Flchtel , 1970). 

Solar Electrons 


and Lin , 1966 ; Lin and Anderson , 1967; 
Anderson, 1969 ; Lin , 1970a, 1971) have been 
Intensively studied. 

The temporal characteristics of relati- 
vistic electrons following the 6 July 1966 
flares are shown in relationship to the 
proton component in Fig. 4. It was found 
that for this flare the electron and proton 
components could be organized in time by 
considerations of particle velocities alone. 
However, it was noted that Lin (1970a) 
showed that the low energy (E > 45 keV) 
electrons apparently arrived first, before 
the more energetic particles. 

An example of a comparison of the elec- 
tron intensity-time profiles for a single 
event Is shown in Fig. 8 ( Lin , 1970b, pri- 
vate communication; Lanzerotti , 1970a). 
These data, from a west limb flare, show 
rather similar time profiles for a wide 
range of electron energies. The similari- 
ties in the temporal profiles are In con- 
trast to those observed in the case of 
protons ( e.g. . Fig. 3; see also Lanzerottl , 
1970 a, for proton temporal profiles mea- 
sured during the same period as the elec- 
tron data of Fig. 8 ). This could be due 
to a more direct propagation of electrons 
to the earth with less interplanetary dif- 
fusion and scattering than In the case of 
protons . 


ELECTRON C0UNTS/S6C 



The first direct observation of solar 
electrons was made from data obtained on 
a balloon flight by Meyer and Vogt ( 1962 ) 
three days after a large flare on 20 July 
1959. They detected highly relativistic 
electrons of energy 100-1000 MeV. Non- 
rela tivistic electrons (E > 45 keV) were 
first measured in interplanetary space by 
Van Allen and Krimigis ( 1968 ) using an in- 
strument flown on Mariner IV. Since that 
time, the time-intensity profiles of rela- 
tivistic solar electrons ( e.g. , Cline and 
McDonald , 1968 ; Simnett et al. , 1959) and 
non-relativistlc solar electrons ( Anderson 
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IV. SPECTRA 

As might be expected, solar flare particle 
spectra have large variations in Intensities 
and spectral shapes both between individual 
events as well as within a single event. 
Representative proton, alpha particle, and 
electron solar flare particle spectra are 
discussed separately. 

Proton Spectra 


The solar proton spectra, particularly 
for higher energies, generally steepen with 
time after the flare. That is, as time 
progresses, relatively few higher energy 
particles as compared to the lower energies 
are present. The proton energy spectra 
measured in several of the large events 
during the last decade were found to fit 
very well a spectral representation with 
a rigidity dependence. This spectral shape 
can be expressed as ( Freier and Webber. 

i » —i 3 


dJ 


dJ. 


dR “ dR ex P[" R /R 0 (t) ] 


( 1 ) 


where R = Mv/c is the proton rigidity. It 
was found that this spectral representation 
was particularly applicable during the de- 
cay phase of an event for protons of ener- 
gies > 20 MeV. 

Six^proton energy spectra obtained during 
several large events of the last solar cycle 
are plotted in Fig. 9. These spectra exhi- 
bit the exponential-in -rigidity spectral 
sha Pe ( Freier and Webber . 1963; . Solar 
cosmic ray spectra such as those of Fig. 9 


PoL,T ery steep compared to the galactj 
cosmic ray spectra (e.g., Fichtfl and 
McDonald . 1967) . v — Gnrel an ° 



Essentially all of the solar flare parti- 
cle spectra taken during the last solar 
cycle were obtained by rocket or balloon- 
based instruments. Frequently, little of 
the time history of an event was obtained. 
Hence, relatively greater emphasis appears 
to have been placed on the spectra of the 
different events. During the present so- 
lar cycle, with essentially continuous 
monitoring of an event* s time profile and 
with the measurement of lower energy par- 
ticles, less emphasis has been placed on 
the individual event spectra. This neglect 
of spectral emphasis partly arises, of 
course, because the spectra changes 
during an event, particularly for the 
events which do not exhibit a diffusive 
temporal profile at the lower energies 
( e.g. , several of the events in Fig. 3). 
Hence, it is impossible to categorize an 
event simply with only one or two spectra. 
However, unlike the more energetic parti- 
cles, and as will be discussed in Section 
V, the decay of low energy particles dur- 
ing events that do have a diffusive char- 
acter appears to be energy independent. 

In this case, a single spectral shape 
would indeed describe much of the event’s 
spectral form. 

A single power-law in energy was fit by 
Lanzerotti (1969c) to the half-hour aver- 
aged proton spectra (E = O.58 to 18.1 MeV) 
measured during the event plotted in Fig. 

6. He found that the spectra became sig- 
nificantly softer during the storm parti- 
cle event on May 24 and for the next two 
and one-half days following the SC on May 
25. At both times the exponent n changed 
from ^1.3 to ~2.0. 

The low energy proton spectra measured 
by Bell Laboratories’ Instruments on 
Explorers 34 and 41 near the intensity 
maximum of several solar flare events in 
the past several years are plotted in Fig. 
10. The spectra are plotted on log-log 
scales to emphasize deviations from simple 
power-law relationships at these energies. 
In particular, the spectrum from the 2 
November 1969 flare has a pronounced peak 
at E ~ 3.5 MeV which may signify particle 
propagation delays from the flare region 
(on the extreme west limb). 

A study of the response of the world- 
wide network of neutron monitors to the 
28 January 1967 event (see Fig. 2) has 
been made by Heristchi and Trottet (1971)* 
They used the global distribution of neu- 
tron monitors as an energy spectrometer to 
determine a possible upper cutoff in the 
energy spectrum of the flare-produced 
protons from this event. They found an 
energy cutoff of 4.3 ± 0.5 GeV. This 
energy is 10^ - 1CP eV lower than that 
predicted by a flare acceleration model 
of Friedman and Hamberger (1969). 


Fig. 9 Solar proton apectra naaaurad In aaveral aolar avanta 
( Fraier and Webber. 1963). 
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o 1730-1800 DAY 145, 1967 UT 

* 1530-1600 DAY 149, 1967 UT ° 0900-1000 DAY 325, 1968 UT 

o 0400-0500 DAY 222,1967 UT • 0330-0400 DAY 103, 1969 UT 

• 1530-1600 DAY 222, 1967 UT « 1600-1630 DAY 306,1969 UT 



Fig. 10 Low energy solar proton spectra measured on Explorers 
34 and 41 near the maximum of several solar events. 


Alpha Particle Spectra 

High energy (E > 30 MeV/nucleon) alpha 
particles were observed in a number of so- 
lar events to have spectra similar to those 
of the protons. Frequently, differential 
rigidity spectra (Eq. l) were applicable 
to both particle species during an event, 
although the e-folding rigidity value R Q 
might at times be different for the two 
species ( e . g . , Biswas et al. , 1963; Dur^a- 
prasad et al M 1968) ^ 

Solar alpha particles were studied over 
a very wide energy range during the 12 
November i 960 solar event. The differen- 
tial alpha fluxes measured between ~31 a nd 
~100 MeV/nucleon by several workers are 
presented in Fig. 11 ( Biswas et al . , 1962; 
Ney and Stein , 1962). Also shown are two 
high energy alpha flux measurements from 
the work of Yates (1964) during the same 
event. * 

The flux measurements of Yates (Fig. 11), 
if expressed on an exponential-in-rigidity 
basis, would fall considerably above what 
would be predicted by an extrapolation of 
the lower energy data ( Ya tes , 1964) . Yates’ 
measurements were challenged by Waddington 
■and Freler ( 1965 ) as perhaps being conta- 
minated by the high fluxes of slow protons 
present during the event. The controversy 
appears still to be unresolved ( Ya tes , 1965 ); 
very high energy measurements of solar al- 
pha particles need to be made during other 
large events . 



Fig. 11 High energy solar alpha particle 

Fluxes measured "by several observers 
following the 12 November i960 solar 
event (from Yates , 1964). 


As noted in Section III, recent years 
have seen an increase in the time resolu- 
tion of low energy solar alpha particle 
observations by satellite. Studies of 
the changes in the alpha spectra during a 
single event have become feasible, although 
little emphasis has been placed on this as- 
pect of the observations. Lanzerotti (1969 
c) studied the power- law exponent of the 
alpha spectra for the May 19&9 event (Fig. 
6). He found that during the period be- 
tween the two sudden commencements (May 24 
and May 25), the alpha particle spectra 
were somewhat harder than that for protons; 
however, after the May 25 SC, both proton 
and alpha- particle power law exponents In 
Eq. (2) were ~2 . The low energy alpha par- 
ticle spectra measured by Bell Laboratories 1 
instruments on Explorer 34 and Explorer 4l 
near the maximum of several solar events of 
the last several years are plotted in Fig. 
12. These spectra, give a representative 
example of low energy alpha spectral shapes 
and intensities. 

Electron Spectra 

Since the number of published observa- 
tions of solar electrons is substantially 
less than for protons, detailed Information 
on electron spectra is less plentiful. Lin 
( 1970 a) has compiled electron spectra from 
four solar events In 1967 as measured by 
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° 1730-1800 DAY 145,1967 UT ° 0900-1000 DAY 325,1968 UT 

* 1530-1600 DAY 149,1 967 UT *0330-0400 DAY 103,1969 UT 



ALPHA PARTICLE ENERGY, MeV/NUCL 


Pig. 12 Low energy solar alpha particle 

spectra measured on Explorers 3 ! and 
£j^near the maximum of several solar 



Fig. 13 Solar electron spectra measured during four events 
In 1967 ( Lin. 1970a}. 


Instruments on Explorers 34 and 35 . These 
spectra are shown In Fig. 13 . Lin finds 
that if he fits the spectra to a power law 
(Eq. 2), the events In Fig. 13 have expo- 
nents n ~ 2 . 3 - 3 . 5 . 

V. PROPAGATION 

After acceleration, the flare-produced 
solar particles must escape from the active 
region and propagate through interplanetary 
space to the earth. Interplanetary space 
Is permeated by the solar magnetic field. 
The nature of this field configuration was 
predicted by Parker (i 960 ) to consist of 
spiral lines emanating from the sun. This 
prediction was subsequently confirmed by 
extensive satellite measurements (e.g., 

Ness et al. , 1964). This spiral inter- 
planetary magnetic field controls much of 
the propagation of the flare particles. 

The guiding center of the solar particles 
tend to follow the spiral nature of the 
field. However, the small-scale irregu- 
larities In the field act as scattering 
centers and perturb, or scatter, the par- 
ticles, moving them to other field lines. 
Extensive theoretical work (not discussed 
here) has been carried out in recent years 
In determining the solar particle diffu- 
sion coefficients due to these random 
scatterings ( e.g. , Jokipli, 1966 , 1967, 

1968; Roelof , 1966, 1968; Hasselmann and 
Wlbberenz , 1968 ; Jokipli and Parker ■ ~1Q5 q) . 


The propagation characteristics of solar 
particles have been reviewed recently 
( Fichtel and McDonald , 1967; Axford , 1970). 
The first considerations of a diffusion 
model for solar particle propagation was 
that of Parker ( 1956 ) and Meyer et al . 
(1956). The broad considerations and 
the development of isotropic diffusion 
theory for solar particles, l.e. , solu- 
tions to a diffusion equation of the form 



5 
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(3) 


have been due to Parker ( 1963 ) . In Eq. 

(3)* n(r,t) is the mean density of solar 
particles with velocity v and k = yv^T = 
3 -Av is the diffusion coefficient and is, 
most generally, a tensor quantity ( Jokipii , 
I 966 ). t is the mean particle "collision 
time" for interaction with the interplanet- 
ary magnetic irregularities and can be 
determined from the interplanetary field 
fluctuations ( Jokipli and Coleman , 1968 ). 
Most commonly, solutions to Eq. (3) have 
assumed k = k (T)r p where T is the parti- 
cle kinetic energy and 0 is time- independ- 
ent . 

Solutions to Eq. ( 3 ) have chiefly con- 
sidered two different boundary conditions. 
The first of these that has been used has 
taken 6=0 and has assumed a perfectly ab- 
sorbing boundary (n=0) at some r = r, > 

1 a.u. The solution to Eq. ( 3 ) at times 
t » r A after the flare yield an expo- 
nential decay for the fluxes 
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(5) 


where the decay time is given as 

2 


T D " 


7T 2 K 


This solution to the model has been util- 
ized by Bryant et al . ( 1962 ) and Hofmann 
and Winckler (19627" in analyzing solar 
particle events. They found that the ab- 
sorbing "boundary" r, was at r ~ 2 a.u. 
Although the 2 a.u. Boundary may indicate 
that the hydromagnetic waves producing the 
interplanetary irregularities are being 
damped out at this distance ( Jokipii and 
Davis , 1969), Axford (1970) has maintained 
that an exponential decay of k with dis- 
tance r will also produce the same results. 

A number of solar events have been fit 
by Krimigis ( 1965 ) using a solution of Eq. 
(3) assuming a radial dependence to the 
diffusion coefficient ( i.e. , P / 0) and no 
boundary r^ ( Parker , 1963 ) * This solution 
can be expressed as 


n(r,t) * f(* o ,P,T) 


h 3/(2-P) n 


exp 


2-p 

r ^ 

(2-e) s 


1 


A plot of ln[n(r,t)t^/^ versus t 

should yield a straight line for the pro- 
per choice of p(< 2 ). Krimigis found that 
for protons in the energy range 5O-5OO MeV, 
good agreement with observations was ob- 
tained for p ~ 1 and A ~ 0.1 a.u. 


Anisotropic Diffusion 

The isotropic solar particle diffusion 
model discussed above is not able to ex- 
plain several important characteristics 
of the solar particles observed at the 
earth. One of these is the direction of 
the non-radial anisotropy of the particles 
measured at the earth. The anisotropy at 
the beginning of events is aligned along 
the spiral field direction, outward from 
the sun ( McCracken , 1963; McCracken et al. , 
1967). (Later in the events the anisotro- 
py becomes much less and the direction 
changes to radial or nearly so (McCracken 
et al. , 1967; Rao et al. , 1969 ) . ) The 
second problem with isotropic models is 
that they can not treat the observations 
that show that the particle fluxes arising 
from flares in the eastern hemisphere of 
the sun tend to increase more slowly to 
maximum intensity than those that origi- 
nate from west hemisphere flares ( e . g . , 
Fichtel and McDonald, 1967; Burlaga , 1967)* 
Reid (1954) has considered a solution to 
the east-west effect by postulating a thin 
diffusing shell around the sun. Particles 
originating from flares in the eastern 
hemisphere would diffuse (isotropically) 
across the solar surface to the spiral 
field lines connecting the sun to the 
earth and then propagate along the field 
lines to the earth (Fig. 14). Reid’s mod- 
el of diffusion across the solar surface 
must be combined with an interplanetary 
propagation model to provide a complete 
description of the particle event as seen 
at the earth. Since inclusion of the 







Fig. lU Illustration depicting solar particle 
diffusion across the solar surface to 
Lhe Interplanetary fin* tube linking 
the sun to the earth ( Reid , 196*0- 


solar-surface diffusion increases the num- 
ber of parameters that can be adjusted, it 
is likely that most diffusive- type obser- 
vations could be fit with such a model. 
Indeed, a solar diffusing layer was one 
of the features included in a recent com- 
putational model for solar flare propaga- 
tion ( Englade , 1971 ). 

The consideration of an anisotropic dif- 
fusion coefficient to solve the east-west 
problem was first made by Axford (1965). 
Burlaga (1967) solved the diffusion equa- 
tion considering particle diffusion trans- 
verse to the spiral interplanetary field 
as well as along it and neglected Reid’s 
diffusion layer around the sun. Expressed 
in spherical coordinates, Burlaga solved 
the equation 





+ 


1 d dn 
2 . 2 q 5cp cp <5p 

r sin 0 ^ 


(7) 


where 1 1 = cos G and Kn, k. , and k are the 
components of the diffusion tensor. 

Burlaga (1967) took the parallel diffu- 
sion coefficient Km to be a constant, in- 
dependent of the radial position, the 
transverse coefficient k. proportional 
to the square of the radial distance 
(k. * r 2 ), and an absorbing boundary at 
rt) 1 > 1 a.u- Solving Eq. (7) with the 
above boundary conditions and assuming n 
to be independent of cp, he obtained quite 
satisfactory fits for a number of differ- 
ent flare events distributed over the so- 
lar disk. The event decay time resulting 
from his solution can be written as 


T 
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Eq. (8) is of the same form as the decay 
time derived for Isotropic diffusion (Eq. 
5) with the isotropic diffusion coeffi- 
cient k q replaced by k . 

Two examples of Burlaga * s fits to parti- 
cle fluxes resulting from flares at two 
separate solar locations are shown in Fig. 
15 . The angle G noted on the figure is 
the angle, measured from the center of the 
sun, between the flare location and the 
location on the sun of the interplanetary 
field line passing through the earth. Fits 
of the model to both the neutron monitor 
observations of the 23 February 1956 event 
and the balloon observations of E > 80 MeV 
protons from the 20 July 1961 event are 
seen to be quite good. 



Low Energy Propagation 

It is clear from data such as those of 
Fig. 3 that at lower energies solar flare 
particles do no t often have diffusive in- 
tensity-time profiles . Substantial modu- 
lation of these low energy particles by 

solar wind discontinuities, shock waves, 
and magnetic field sector boundaries must 
be occurring. Although there are events 
where the low energy particles exhibit 
diffusive-type profiles, the applicability 
of an anisotropic diffusion model such as 
Burlaga's to these observations is highly 
suspect ( Forman , 1970 ). Forman has main- 
tained that the "equilibrium" anisotropy 
present during the decay phase of an event 


( McCracken et al. , 1967) and the evidence 
that the diffusion coefficient becomes 
small at low energies ( Jokipii and Coleman . 
1968) indicate that solar wind convection, 
and the resulting particle energy loss, is 
an important mode of low energy particle 
propagation. (The anisotropic diffusion 
model of Burlaga (as well as the isotropic 
models) considers convection effects to be 
negligible for the higher energy particles, 
and rightly so.) Forman (1971) cites as 
further evidence for the importance of 
convection the fact that the reported de- 
cay times for both protons and alphas were 
essentially energy-independent for 1-20 
MeV/nucleon particles after the 28 May 
1967 flare (Lanzerotti, 1969a). 

Forman (I 9 ?l) has solved the Fokker- 
Planck equation first derived by Parker 
(1965) for particle transport including 
convection and diffusion: 

^ i >*>) - -4 

■ 3 £ (9) 

Here V is the solar wind velocity, T is 
the particle kinetic energy, and a = 

(T+ 2 Mc )/(T+Mc ). The diffusion models 
discussed above all neglected the terms 
in Eq. (9) containing the solar wind velo- 
city V. Forman obtained an analytic solu- 
tion to Eq. (9) assuming that k. = K]_rS 
Km = K^r , and that there was a diffusing 
boundary at r - r, . Formats model pre- 
dicts very well the equilibrium residual 
anisotropy during the decay phase of an 
event as observed by McCracken et al. 

(1967) as well as the - magnitude (~ 14-18 
hours) of the energy-independent decay 
time for both alphas and protons as ob- 
served by Lanzerotti (1969a). 

From her solution to Eq. (9) Forman ( 1971 ) 
has shown that in the decay phase of the 
event 
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where the decay time is given as 
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In Eq. (11) 


is 


Bessel functiSh^of order T|. 


found that for r^ = 2.3 


the first zero of the 
Forman has 


a ,u . 


(a representa- 


tive value determined from the model fits 
of - U 1 ? y 3 h 1967 b r D has a broa d maximum 

aLt^Tn?0 h °2 rs £?Vll = between - 3 - 10 19 
and 3 10 ^ u cm c sec 1 (r^ason&ble values for 

| aa determined from the power spectra of 
tne interplanetary magnetic field near the 
earth by Jokipii and Coleman . 1968) For- 
k’s model has recently been applied 
cessfuiiy to the low energy pro?on oSser- 
vations from ^ the 7 June 1969 event ( Murray 
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The energy-independence of the decay 
times for both protons and alpha particles 
during a diffusive-type event is shown in 
Fig. 16 for the 13 April 1969 event (Lan- 
zerotti and Graedel, 1970). The intensity- 
time "prof iles of the proton fluxes in the 
O .56 < E < 0.60 MeV channel is shown as an 
insert in the figure. This event, probably 
originating from a flare behind the east 
limb, demonstrated a diffusive-type appear- 
ance even in the lowest energy-channel mea- 
sured. This was quite unlike the 28 May 
1967 diffusive event where energetic storm 
particles greatly enhanced the lower energy 
proton fluxes ( Lanzerottl , 1969 a). Also 
plotted In Fig. 16 are”the decay times for 
E > 10, >30, and >60 MeV protons measured 
by the solar particle monitoring experiment 
on the same satellite ( Solar Geophysical 
Data, 1969). The decay time varied from 
~28 hours at O .58 MeV to ^19 hours. at 60 
MeV. Over the range O .58 MeV to 20 MeV, 
the decay time decreased by only ~4 hours. 



Pig. 16 Proton and alpha particle decay times following the 

13 April 1969 Bolar event. The intensity -time profile 
for the O.56 £ E £ 0.60 MeV proton channel is shown 
as an inset to the figure ( Lanzerottl and Graedel, 1970). 


VI. FLARE PARTICLE EFFECTS 

Two consequences of solar flare particle 
effects are discussed below. The first of 
these is the effect of energetic flare par- 
ticles In producing polar-cap cosmic noise 
absorption and the detection of this en- 
hanced absorption by riometer techniques. 
The second is the effect of solar parti- 
cles, penetrating into the outer magneto- 
sphere, on satellite solar cell lifetimes. 


Riometer Absorption 

The first indication of the production of 
enhanced ionosphere ionization by solar 
flare particles was the strong absorption 
of cosmic radio noise in the polar cap re- 
gions that Bailey (1957) correlated with 
the flare of 23 February 1956. Since that 
time, enhanced riometer absorptions in the 
auroral and polar cap regions during solar 
events have been studied as basic geophysi- 
cal phenomena and as diagnostic tools for 
studying soldr and magnetospheric processes 
(e . g. , Bailey , 1964; Reid , 1970). Indeed, 
the significance of energetic storm parti- 
cles was first outlined by Axford and Reid 
(1963) using riometer data. 

Through the work of Potemra et al. (1967. 
1969 , 1970 ) good agreement has been achieved 
in calculating the expected riometer re- 
sponse from a measured incident solar flux. 
Potemra and his collaborators have calcu- 
lated the expected total absorption A at a 
radio wave angular frequency cjd from the 

formula 


A(dB) = 


I.l6xl0 6 is 
m 



(12) 


obtained from the theory of Sen and Wyller 
(i 960 ). In Eq. (12) ovr is the angular 
gyro frequency, dh is The increment of 
ionization height in 10-km units, C^/o 
an integral function, n is the electron 
density and v is the mean electron colli- 
sion frequency. Using specific ionization 
rates due to G. W. Adams and Adams and 
Mas ley ( 1965 ) and v and recombination 
coefficients deduced from the September 
1966 event, P otemra et al. (1970) have 
predicted the observed absorption for high 
latitude riometer observations during a 
number of 1967 PCA events. Their calcula- 
tions, using satellite measurements of the 
solar proton fluxes over the polar caps, 
are compared in Fig. 17 to the observed 
riometer day and night absorption measure- 
ments following the 28 January 1967 solar 
event. The agreement is quite good for 
both the day and the night observations, 
(it is interesting to compare these lower 
energy proton observations of Fig. 17 with 
the neutron monitor profile for the same 
event in Fig. 2 . ) 

Several authors ( Van Allen et al. , 1964; 
Juday and Adams , 1969; Reid ", 1969. 1970) 
have used the empirical relation 

(F)^ = RxA (13) 


to relate the Integral fluxes F of protons 
above some energy E . to the riometer ab- 
sorption A at a givenirequency . In Eq. 
( 13 ). R is a constant, dependent only upon 
E . . Potemra and Lanzerottl (1971). using 
s<5l§r proton data from the synchronous 
equatorial ATS-1 satellite, deduced R as 
a function of E , from the 30 MHz riometer 
absorption obse?v§d at Byrd during the 28 
January 1967 event. (Byrd (le-7) is at 
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Solar Cell Damage 



Pig. 17 Comparison of observed day and night rlometer absorption 
and the calculated absorptions using polar cap-average 
solar proton data from satellite 1963 - 38 C during the 
28 January I 967 event ( Potemra et. al. . 1970). 


nearly the same latitude as ATS-1 (Ia- 6.4) 
but three hours earlier In local time.) 
They found Eq. (13) to be an excellent 
fit to the data, essentially Independent 
of the value of E , . The values of 
R(E miJ the y found for E , between 5 
and 5“ MeV are plotted iffi 1 j?lg. 18 as a 
function of E , . Such R-values should 
be quite useful for the new solar proton 
event classification scheme (Shea and 
Smart , 1970). 


Solar protons appear to have ready access 
to the outer regions of the magnetosphere 
other than through the polar cap regions. 
The solar particles at synchronous altitude 
are observed to have essentially the same 
intensities and spectra as the particles 
In Interplanetary space for protons as low 
as 1 MeV in energy ( Lanzerottl , 1968, 1970b; 
Paulikas and Blake , 1969 ) . These low ener- 
gy protons could cause significant damage 

to unshielded solar cells on a synchronous 
satellite. It was pointed out by Lanzer- 
otti (1969b) that the damage to unshielded 
cells from relatively low Intensity solar 
events could dominate the normal synchron- 
ous altitude radiation (predominantly 
electrons) in producing damage. Indeed, 
anomalous, step-like changes in the short 
circuit current of unshielded cells in the 
ATS-1 solar cell damage experiment are ob- 
served in conjunction with solar flare 
events ( Waddel , 1968). 

An example of the solar proton damage to 
an unshielded solar cell on ATS-1 during 
the May 1967 events (Fig. 6 ) is shown in 
Fig. 19. The two lower bar graphs show 
the values of the short circuit current 
(in ma) measured each day for two 10 O-cm 
n-on-p type solar cells (R. C. Waddel, 
private communication). One cell was un- 
shielded while the other had a 1 mil shield 
of 7740 glass. At the top are plotted the 
daily average of the integral half-hour 
average proton fluxes (E > 2.4 MeV) mea- 
sured by the Bell Laboratories experiment 
on ATS-1. (No data were received on davs 
143, 144, and 147.) 


E 8 I— - 



. lfi Factor R. Empirically relating the E , MeV 

rlometer absorption to the square root MIN 
of the Integral proton fluxes plotted 
as a function of the lover 

energy limit on the Integral fluxes. 

R vas determined fro. the ?P January 
1W event ( Fount™ • 

1971 1 . 


SHORT CIRCUIT CURRENT ,ma. LOG FLUX! >2. 4 MeV) 
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Prior to the major interplanetary en- 
hancement on day 145, the proton fluxes 
at ATS- 1 were very low. The day after 
the large ATS-1 enhancement on day 145 
the short circuit current in the unshield- 
ed cell decreased sharply whereas only a 
slight decrease in the current was observed 
in the shielded cell. 

Since a 1 mil shield will stop ~0.5 MeV 
protons, the data of Fig. 19 are indicative 
that protons with energies as low as this 
were producing the most significant damage. 
This could be due both to the penetration 
of 0.5 MeV solar protons to the synchron- 
ous orbit as well as to the fact that the 
magnetosphere boundary was pushed within 
the ATS-i orbit for periods of time on 
days 145 and 146. 

Even if shields were provided for syn- 
chronous satellite solar cells, the fluxes 
of low energy (0.5-3 MeV) solar protons 
penetrating to this altitude could still 
play an important role in the damage con- 
siderations. This is because in the manu- 
facturing process the shields are often 
not deposited uniformly over the cell 
surfaces, leaving a small fraction of 
some cells uncovered (R. C. Waddel, pri- 
vate communication). Such partially un- 
shielded cells would then be subjected 
to unexpected damage by the low energy 
protons . 
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